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O^. Abstract. In this geometrical approach to gravitational lensing theory, we apply the 

^ ' Gauss-Bonnet theorem to the optical metric of a lens, modelled as a static, spherically 

^^ ' symmetric, perfect non-relativistic fluid, in the weak deflection limit. We find that the 

focusing of the light rays emerges here as a topological effect, and we introduce a new 
method to calculate the deflection angle from the Gaussian curvature of the optical 
metric. As examples, the Schwarzschild lens, the Plummer sphere and the singular 
l/~j I isothermal sphere are discussed within this framework. 
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1. Introduction 

The deflection of liglit by gravitational fields has been studied with great interest 
in astrophysics as well as in theoretical physics. Fundamental properties such as 
Fermat's principle for Lorentzian manifolds, conditions on image multiplicity and 
caustics in spacetime have been discussed in a fully relativistic setting (e.g., see P and 
references therein). In the astrophysical context, however, an impulse approximation 
with piecewise straight light rays in fiat space has proven useful since deflection angles on 
cosmological scales are very small (for a comprehensive introduction, see e.g. [2] or [3] pp 
272-95 and references therein). Despite their different premisses, both treatments have 
yielded mathematically interesting, general properties which depend on topology. In 
particular, image counting theorems like the odd number theorem have been established 
with different versions of Morse theory both in the spacetime lensing and impulse 
approximation frameworks [1]. 

In this article, we would like to present another approach to gravitational lensing 
theory which emphasizes global properties. Specifically, we consider the astrophysically 
relevant weak deflection limit not in the impulse approximation but treat light rays 
as spatial geodesies of the optical metric, and use the Gauss-Bonnet theorem. This 
approach has previously been applied to lensing by cosmic strings [5], and we extend 
it here to static, spherically symmetric bodies of a perfect fluid as simple models for 
galaxies acting as gravitational lenses. It turns out that the focusing of light rays is, from 
this point of view, essentially a topological effect. Hence we find, rather surprisingly, 
that the deflection angle can be calculated by integrating the Gaussian curvature of the 
optical metric outwards from the light ray, in contrast to the usual description in terms 
of the mass enclosed within the impact parameter of the light ray. To illustrate, we 
discuss how this works for three well-known models, namely the Schwarzschild lens, the 
Plummer sphere and the singular isothermal sphere. 

The structure of this article is therefore as follows. In section 2 we give a brief 
review of the Gauss-Bonnet theorem and introduce two constructions which will be 
used to investigate the lensing geometry. The optical metric and its Gaussian curvature 
for static, spherically symmetric systems of a perfect fluid is discussed in section 3, 
followed by the application to the three lens models mentioned above in section 4. 

With regard to conventions, we use metric signature (—,+,+,+), Latin and Greek 
indices for space and spacetime coordinates, respectively, and set the speed of hght 
c = 1. G denotes the gravitational constant as usual. 

2. Gauss-Bonnet theorem and lensing geometry 

The Gauss-Bonnet theorem connects the intrinsic differential geometry of a surface with 
its topology, and this will be the main tool in our exposition. Let the domain {D, x, g) 
be a subset of a compact, oriented surface, with Euler characteristic x and a Riemannian 
metric g giving rise to a Gaussian curvature K. Furthermore, let dD : {t} ^ D he 
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Figure 1. Weak deflection lensing geometry. Two geodesies 71 and 72 from the 
source S to the observer O are deflected by a lens with centre at L. Di and D2 are 
two domains with boundary curves 7^ and 7p as discussed in the text. 

its piecewise smooth boundary with geodesic curvature k, where defined, and exterior 
angle ai at the ith vertex, traversed in the positive sense. Then the local and global 
versions of the Gauss-Bonnet theorem (see, e.g., [6] pp 139, 143) can be combined to 
give 

KdS+ f Kdt + y2ai = 2nx{D). (1) 

'D JdD ■ 

Next, consider a smooth curve 7 : {t} ^ D of unit speed such that ^'(7,7) = 1, and let 
7 be the unit acceleration vector. This vector, then, is perpendicular to 7 and hence 
spans a Frenet frame together with 7. The geodesic curvature of 7 is therefore given by 
(see, e.g., [6] p 138) 

fi: = (7(V^7,7), (2) 

which clearly vanishes if, and only if, 7 is a geodesic. We will now describe two specific 
setups for the domain, Di and D2, which will represent the weak defiection lensing 
geometries discussed in section 4, and are shown in figure [H 

Firstly, let Di be bounded by two geodesies 71,72 intersecting in two vertices, the 
first being the source S, and the second the observer O. Hence ^(71) = ^(72) = except 
at the source and observer where the interior angles are 63 = ii — as and do = t^ — oiq, 
respectively. Di contains the lens centre L, and both S and O are assumed to be very 
distant from L such that 9s and Oq are small and positive. If the lens centre is non- 
singular, then Di is simply connected so that x{Di) = 1- On the other hand, if L is 
singular, then x{Di) = and there is a contribution to dDi given by the boundary 
curve 72, enclosing L. Hence the Gauss-Bonnet theorem ([1]) yields 

^s + 0o= [ [ KdS (3) 

in the non-singular case, and 

0g + 0^ = 2n+ (f Kdt+ f f KdS (4) 
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in the singular case. 

Secondly, let D2 be a simply connected domain such that L ^ D2 and Di nD2 = 71. 
Furthermore, D2 be asymptotically flat such that S and O are in this regime. Thus 
dD2 consists of a geodesic 71 from S" to O, and the perimeter curve ■jp of a circular 
sector centred on L intersecting 71 in 5* and O at right angles. If (p is the angular 
coordinate centred on L, then 7^ has angular range tc + 6 where 6 is the small and 
positive, asymptotic deflection angle of 71. Therefore we can use ([1]) to write 

^ «:(7p)^d(/)-7r = -// KdS, (5) 

again using ^(71) = 0. By construction, the calculation of the deflection angle according 
to ([5]) is obviously independent of any singularity that may occur at L in Di. 

3. Optical metric and its Gaussian curvature 

We shall now consider null geodesies deflected by a static, spherically symmetric massive 
body consisting of a perfect fluid. In the present context, this is to be thought of as the 
stellar fluid of a galaxy acting as a gravitational lens. The corresponding line element 
written in Schwarzschild coordinates x'^ = (t,r,{},(f)) is given by (e.g., see [3j p 375) 

ds^ = g^u(ix'^dx'^ 

= - exp[2A{r)]dt^ + exp[2B{r)]dr^ + r^ {d9^ + sin^ M0^) . (6) 

By spherical symmetry, we can assume, without loss of generality, that the null geodesies 
with ds^ = are in the equatorial plane -^ = n/2. All images are therefore collinear 
with the lens centre L in this case. Spatial projection of the null geodesies yields the 
light rays, and these are geodesies of the optical metric g!^^ = gmn/{—goo), by Fermat's 
principle. The geometry of the optical metric, also known as optical reference geometry, 
is useful for the study of gravitational and inertial forces in general relativity |7] . It is 
convenient to introduce a radial Regge- Wheeler tortoise coordinate r* defined by 

dr* = drexp[i?(r) — v4(r)], 

so that the line element of the optical metric becomes 

dt^ = C^dx'^dx'^ = exp[2(B - A)]dr^ + exp{-2Ayd(f)^ 

= dr*2 + f{r*yd(f)^ (7) 

from ([6]), where /(r*) = exp[— A(r)]r, r = r{r*). It is clear, then, that the equatorial 
plane in the optical metric is a surface of revolution when it is embedded, where possible, 
in M^. Its (intrinsic) Gaussian curvature K (e.g., see [6J p 66) can be found from ([71) 
and expressed in terms of the Schwarzschild radial coordinate r, 

1 dV 



K 



ZdAdB IdA IdB d'^A' 
\ dr dr r dr r dr dr^ ' 



Applications of the Gauss-Bonnet theorem to gravitational lensing 5 

One can now specialize to the perfect fluid case. Then the energy-momentum tensor is 
rptiu _ cliag(p,p, p,p) in the local flat metric, where p = p{r) denotes the density and 
p = p{r) the pressure of the lens model. Now from Einstein's field equations (e.g., see 
[31 PP 376-7), 

exp[-25(r)] = 1 - -^ where p{r) = AnG / p{r'y^dr'. (9) 

^ Jo 

The conservation of the energy-momentum tensor, together with the field equations, 

gives rise to the Tolman-Oppenheimer-Volkoff equation 

dp (p + p) {p + AnGr^p) 



dr r2 (l-f) 

which expresses the hydrostatic equilibrium, and we also have 
dA 1 //i 



(10) 



r 

Hence using ([8]), the Gaussian curvature is 

2pexp[2{A-B)] (^ 3fi AnGr^ 
K = ■ —-^ — < 1 



+ 4:7TGpr). (11) 



2p\ ^ f^ 3p ^_^ 2 



Since we are presently interested in weak deflections only, one can limit the discussion 
to a non-relativistic stellar fluid subject to the collisionless Boltzmann equation (e.g. 
^). But this means that the pressure may be neglected in ([H]). To see this, recall that 
non-relativistic kinetic theory, for instance in the case of an isotropic velocity dispersion 
0"^, implies that p = po"^/3 but a^ <c; 1. Now the area element of the equatorial plane 
in the optical metric is 



dS" = ^/\detg°P^\drd(f) = exp[5(r) - 2A{r)]rdrd(f) 



from ([7]). Together with (jH]), the integrand of the Gaussian curvature therefore reduces 
to 



KdS = '^ 



r 



(1-^) 



2/x\3/2 



3/i AnGr^p / 2p 
2r p \ r 



drd0. (12) 



Notice that for reahstic lens models, where p -^ po > as r — > 0, and p —>■ poo < cx) as 
r — > cxD, the Gaussian curvature must change sign at some radius in the equatorial plane. 
For if r ^ oo, then p has to decrease faster than r~^ for p to be asymptotically finite, 
so -ft" < follows immediately from flT2l) . Conversely, if r — ;> 0, then p{r) -^ ATiGp^r^ /?> 
by dHD, so the term in the square bracket in (I12p tends to —2, and the result follows. 
In the next section, we shall discuss three particular lens models as examples. 
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4. Applications to lensing models 

4.I. Schwarzschild lens 

The Schwarzschild lens can be characterized by fi{r) = const, and p{r) =0, r > so 
that outside the event horizon at r = 2/i, 

KdS = ^-^^-^ f 1 - ^"j < (13) 

from flT2]) . The optical Schwarzschild metric has therefore negative Gaussian curvature, 
and an embedding diagram illustrating this can be found e.g. in f?]. Now, in light 
of this fact, it might appear surprising at first glance that light rays can be focused. 
For consider a geodesic ■jit) separated from a neighbouring geodesic by a Jacobi field 
Y{t) = y{t)j, then the geodesic deviation equation becomes (see, e.g., [6] p 102) 

|| + ^(7)l/ = 0, 

implying that the light rays should diverge. However, given a domain Di as described 
in section 2, the focusing of two light rays at the observer's vertex O is made possible 
by the non-trival topology of Di such that x(-Di) = because of the event horizon. It 
is convenient to use the circular photon orbit at r = 3/i as inner boundary curve 7^ 
because its geodesic curvature vanishes. Then equation (jl]) applies, 

9s + 9o = 27i+ f f KdS, 

which can clearly be fulfilled for small, positive 9s, 9o despite the negative Gaussian 
curvature. The topological contribution is therefore essential to the focusing, and it is 
in this sense that gravitational lensing can be understood as a topological effect. In 
fact, its global nature is also borne out by equation ([5]). Here, the deflection angle is 
calculated by integrating the Gaussian curvature in the domain D2 outwards from the 
light ray, as opposed to the usual treatment where the deflection is determined by the 
mass enclosed within the impact parameter (e.g. [2] p 231 and [3] p 277). Now since the 
optical Schwarzschild metric is asymptotically Euclidean, we can take K{'yp)dt/d(f) = 1 
on the circular boundary of D2. Moreover, because of the weak deflection limit, we 
may assume that the light ray is given by r{t) = 6/sin0 at zeroth order with impact 
parameter b ^ 2/i. Hence (|5]) together with flT3l) implies 

6 = - f f KdS -II ^drd0 

J Jd2 Jo Jb/sin(t> ^ 

- 1 (-) 

This is the well-known formula for the Schwarzschild deflection angle in the weak limit 
(e.g., [2] p 25), as required. 

Finally, it should be remarked that the negative Gaussian curvature of the optical 
metric is a rather general feature of black hole metrics, as is the closed photon orbit 
forming part of 9Z?2- The existence of photon spheres for static, spherically symmetric 
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metrics was first discussed in [H], and more generally, in terms of an energy condition, in 
|10j . A fuller discussion of the Schwarzschild case, including applications of the Gauss- 
Bonnet theorem to geodesic triangles, can be found in the forthcoming article [TT] . 



4-2. Plummer model 

The next example is one of the simplest realistic (in the sense of section 3) models for 
a stellar system, the Plummer sphere [12]. It is a polytrope with density distribution 
(e.g., iSj p 225) 

21 -5/2 



p{r) 



Po 



1 



(15) 



where po is the central density and ro a scale radius, so that the mass parameter defined 
in (19]) becomes 



fi{r) 



r 
Poc I — 



r 
1+1 - 

^0 



-3/2 



(16) 



The corresponding metric is therefore an interior solution with finite /ioo proportional to 
the total Newtonian mass of this model in the non-relativistic limit where /ioo/''"o ^ 1- 
We can now compute the Gaussian curvature integrand from flT2l) to first order in fioo/ro, 
to obtain 

2/ioo M /^ 3 



KdS 



dud4> + O ( ^ 



(17) 



ro (1+^2)3/2 y 1+u'y V ,0 

where a dimensionless radius u = r/r^ has been introduced for convenience. Since fi is 
finite at large radii, the optical metric of the Plummer sphere approximates the optical 
metric of the Schwarzschild lens and its negative Gaussian curvature at large radii. Of 
course, there is no event horizon in this model, so the lens centre L is non-singular, 
and the domain Di is simply connected. Then the relevant form of the Gauss-Bonnet 
theorem ([3]) requires that K > somewhere in Di for focusing to be possible. Indeed, it 
can be seen from (fT7|) that the leading term of the Gaussian curvature changes sign at 
u = y2, and this is in agreement with the more general conclusion mentioned in section 
3. 

We can now calculate the deflection angle for the Plummer model by considering 
the domain D2 as in the previous example. From ([S]) and flTTl) . 




KdS 



D2 



2/Xc 



4/ioo h 



1 + 



fe/(rosin0) 
-1 



'^O 



u 



3u 



(1+^2)3/2 (1+^2)5/2 



dud(j) 



after suitable substitutions. This is the weak deflection angle as expected (e.g., [2] p 
245), and we recover the Schwarzschild deflection angle (IT^ in the appropriate limit 
b > ro. 
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4-3. Singular isothermal sphere 

The last model discussed here is the singular isothermal sphere. Fully relativistic, 
isothermal solutions of the Tolman-Oppenheimer-Volkoff equation exist [I3], and the 
density distribution of the singular isothermal sphere is, in fact, similar to the global 
monopole [H]. However, since we are presently concerned with the weak deflection limit, 
it suffices to consider the non-relativistic, singular isothermal sphere with the following 
density distribution (see, e.g., |3| pp 289-90 and [8] p 226-8) and mass parameter, 

,2 



P(^) = o-^^2 ^ ^('') = 2cr r. 



a 
2^r2 

using ([9]), where cr^ ^ 1 is again the isotropic velocity dispersion of the stellar fluid. 
This model is a solution of the Tolman-Oppenheimer-Volkoff equation (fTOl) in the non- 
relativistic limit, that is, to first order in o"^. Then ( IT2l) shows that the Gaussian 
curvature vanishes for r > in this limit, so that the equatorial plane in the optical 
metric is a cone with a singular vertex at r = 0. Another way to see this is to 
isometrically embed the optical metric ([7]) in M.^ with cylindrical coordinates {z, R, 0) 
such that 

dt^ = exp[2(5 - A)]dr^ + exp{~2A)r^d(f)'^ 
= dz{ry + di?(r)2 + R{rfd(f)\ 

Using the equations ([9]), (ITT|) and (IT8|) . one can find the coefficients of the spacetime 
metric (EI), 

exp(2v4) = C'-V^"'/(i-^"'), 

exp(25) = (1 - 4:a^y\ 
where C is a non-zero constant. Then the embedding yields 

R{r) = Ct^^-^^'^/^'-'^'\ 
and the equatorial plane in the optical metric is described by a cone in R'^ of the form 

9\ 1/2 
1 _ 9£!. 



^ = ^8^2^- — ^A^R 

1 — 6(7^ 

in the non-relativistic limit. This corresponds to a deficit angle of A ^ %'Ka'^ . 

Now in order to find the defiection angle, we can again apply Gauss-Bonnet in the 
form of equation (j5]) which reads 

r"^^ dt 

/ fi:(7p)— d0-7r = O (19) 

since K = Q in D2- The Gaussian curvature of the circular perimeter curve 7p can be 
calculated directly from ([2]), 



K{-ip)dt= {l-Aa^f'^(l 



dA\ 
dr J 
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using (II ip . The leading term of the deflection angle is therefore 

from (TT9l) . in agreement with the standard treatment (e.g., J3J p 290). This shows that, 
from this point of view, the constant deflection angle of the singular isothermal sphere 
comes essentially from the deficit angle of the conical optical metric, similar to cosmic 
string lensing (e.g., see [5]). 

5. Concluding remarks 

In this article, we have introduced a geometrical approach to gravitational lensing theory 
different from the spacetime and impulse approximation treatments. By applying the 
Gauss-Bonnet theorem to the optical metric, whose geodesies are the spatial light rays, 
we found that the focusing of light rays can be regarded as a topological effect. We also 
gave a new expression ([5]) to calculate the deflection angle by integrating the Gaussian 
curvature of the optical metric outwards from the light ray, in the weak deflection limit. 
The lens models considered were given by static, spherically symmetric bodies of a 
non-relativistic, perfect fluid, and we discussed as examples the Schwarzschild lens, the 
Plummer sphere and the singular isothermal sphere. 

It would therefore be interesting to see whether this approach could also be extended 
and be fruitfully applied to lenses without spherical symmetry, where images are no 
longer coUinear with the lens centre, or to the relativistic strong deflection limit. 

Acknowledgments 

MCW would like to thank Claude Warnick for useful discussions, and the Science and 
Technology Facilities Council, UK, for funding. 

References 

[1] Perlick V 2000 Ray Optics, Fermat's Principle, and Applications to General Relativity (Berlin: 

Springer- Verlag) 
[2] Schneider P, Ehlers J and Falco EE 1999 Gravitational Lenses (Berlin: Springer- Verlag) 
[3] Straumann N 2004 Genera? i^eZoiiuiij/; With Applications to Astrophysics {Qerlui: Springer- Verlag) 
[4] McKenzie R H 1985 A gravitational lens produces an odd number of images J. Math. Phys. 26 

1592-96 
Fetters A O 1995 Multiplane gravitational lensing. I. Morse theory and image counting J. Math. 

Phys. 36 4263-75 
[5] Gibbons G W 1993 No glory in cosmic string theory Phys. Lett. B 308 237-39 
[6] Klingenberg W 1978 A Course in Differential Geometry (New York: Springer- Verlag) 
[7] Abramowicz M A, Carter B and Lasota J P 1988 Optical reference geometry for stationary and 

static dynamics Gen. Rel. Grav. 20 1173-83 
[8] Binney J and Tremaine S 1987 Galactic Dynamics (Princeton: Princeton University Press) 
[9] Atkinson R d'E 1965 On light tracks near a very massive star Astron. J. 70 517-23 
[10] Claudel C-M, Virbhadra K S and EUis G F R 2001 The geometry of photon surfaces J. Math. 

Phys. 42 818-38 



Applications of the Gauss-Bonnet theorem to gravitational lensing 10 

[11] Gibbons G W and Warnick C M 2008 Hyperbolic space and no hair theorems {forthcoming) 
[12] Plummcr H C 1911 On the problem of distribution in globular star clusters Mon. Not. R. Astron. 

Soc. 71 460-70 
[13] Bisnovatyi-Kogan G S and Zcl'dovich Ya B 1969 Models of clusters of point masses with large 

central redshifts Astrophys. 5 105-9 
[14] Barriola B and Vilenkin A 1989 Gravitational field of a global monopolc Phys. Rev. Lett. 63 341-43 



